Directed polymer in a random medium of dimension 1 + 3 : 
multifractal properties at the localization/delocalization transition 
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We consider the model of the directed polymer in a random medium of dimension 1 + 3, and 
investigate its multifractal properties at the localization/delocalization transition. In close analogy 
with models of the quantum Anderson localization transition, where the multifractality of critical 
wavefunctions is well established, we analyse the statistics of the position weights WL(r) of the end- 
point of the polymer of length L via the moments Y q (L) = ^2jt[wL,(r)] 9 . We measure the generalized 



exponents r(q) and f(q) governing the decay of the typical values Y q typ (L) = e llxY « {L) ~ L~ T{q) and 
disorder-averaged values Y q {L) ~ L~ T ^ respectively. To understand the difference between these 
exponents, r(q) 7^ f (q) above some threshold q > q c ~ 2, we compute the probability distributions of 
y = Y q (L)/Y q typ (L) over the samples : we find that these distributions becomes scale invariant with 
a power- law tail 1/y 5 . These results thus correspond to the Ever-Mirlin scenario [Phys. Rev. 
Lett. 84 , 3690 (2000)] for the statistics of Inverse Participation Ratios at the Anderson localization 
transitions. Finally, the finite-size scaling analysis in the critical region yields the correlation length 
exponent v ~ 2. 



I. INTRODUCTION 



*^ \ The notion of multifractals is now important in various areas of physics (see for instance 0, 0, H, 0, H, S 
£^ ■ and references therein). For classical systems with frozen disorder of interest here, the idea that multifractality is 
present at criticality has been mostly investigated for correlation functions in two-dimensional diluted ferromagnets 
H, [t| Qil [HI 03) spin-glasses and random field spin systems [H, Qj], 0j|. For disordered quantum spin-chains, the 
statistics of critical correlation functions is described by "multiscaling" , which is even stronger than multifractality 
[l6| . For quantum localization models, the multifractality of the critical wavefunction shows up through the statistics 
of inverse participation ratios (I.P.R.s) [Tt], EH- Many results are now available for the Anderson localization transition 
in d = 3 [l^, T2TI . |22| , the integer quantum Hall transition [23|, [24| Dirac fermions in a random magnetic field [2|| , 



and power-law random banded matrices [26[ . Connections have been also established with the scaling properties of 
the correlation functions [27| and with the time evolution of wave packets 28] . More recently, it was realized that 
typical and disorder-averaged I.P.R.s can actually lead to two different multifractal spectra as a consequence of the 
t~^. ' broadness of their probability distributions [2l|, |2g, l29j . 

In this paper, we consider the localization/delocalization transition of the directed polymer in a random medium 
of dimension 1 + 3 (see the review (3fjj and references therein) , to investigate whether some multifractality in present 
^3 i at criticality, in analogy with the quantum localization models quoted above. 

The paper is organized as follows. In Section [TTJ we introduce the directed polymer model and the observables 
displaying multifractal behavior at criticality. We then describe our numerical results concerning the generalized 
dimensions D{q) and D(q) that govern typical and averaged values (Section |TTT]), the singularity spectrum (Section 
lIVp . and the probability distributions over the samples ( Section |Vj). Finally, we present in Section PVTl the finite-size 
Q . scaling analysis in the critical region. Section lyill contains our conclusions. 



II. MODELS AND OBSERVABLES 
A. Model definition 

The random bond directed polymer model [3fj| is defined by the following recursion relation for the partition function 
on the cubic lattice in d — 3 

2d 

Z t+1 (r) = J2 e-M** ^ Z t (r + e, ) (1) 

3=1 
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The bond energies et(r + are random independent variables drawn from the Gaussian distribution 



In this paper, we consider the following boundary conditions. The first monomer is fixed at r = 0, i.e. the initial 
condition of the recurrence of Eq. fl} reads 

Zt=o(f) = S pg (3) 

The last monomer is free, i.e. the full partition function of the polymer of length L is then obtained by summing over 
all possible positions r at t = L 

Zl ot = Y,Z L (r) (4) 

r 

The phase diagram of this directed polymer model as a function of space dimension d is the following 30]. In 
dimension d < 2, there is no free phase, i.e. any initial disorder drives the polymer into a strong disorder phase, 
where the order parameter is an 'overlap' [3~il . [33 . [H, [U, In dimension, d > 2, there exists a pha se transition 
between the low temperature disorder dominated phase and a free phase at high temperature [H, I37ll38l| . This phase 
transition has been studied exactly on a Cayley tree [3l| . In finite dimensions, bounds on the critical temperature T c 
have been derived [13, IH, H(| : T (d) < T c < T 2 (d). The upper bound T 2 (d) corresponds to the temperature above 
which the ratio Z\j (Z^) 2 remains finite as L — > oo. The lower bound To corresponds to the temperature below which 
the annealed entropy becomes negative. On the Cayley tree, the critical temperature T c coincides with To [3l|. In 
finite dimensions however, we have argued in [4l| that T c coincides with T 2 , and our recent numerical simulations [42| 
are in agreement with the numerical value given in [39| for T 2 (d — 3) : 

T c = T 2 (d = 3) = 0.790.. (5) 

The numerical results given below have been obtained using polymers of various lengths L, with corresponding 
numbers n s (L) of disordered samples with 

L = 6,12,18,24,36,48,60,72,84,96 (6) 
n s (L) = 10 s , 10 7 , 2.10 6 , 8.10 5 , 2.10 5 , 5.10 4 , 3.10 4 , 2.10 4 4.10 4 , 2.10 4 (7) 

In the following, A denotes the average of A over the disorder samples. We also define (AA) 2 = A 2 — A" . 

B. Notion of multifractal statistics at criticality 

In this paper, we will focus on the statistical properties of the weights 



Zl(t) 

w L yr) = 

normalized to (Eq. 0|) 



= (8) 



zi 



X>L(f) = l (9) 

r 

and study whether they present some multifractal properties at criticality. We thus consider the following moments 
of arbitrary order q 

F,(L) = X>!(r) (10) 

r 

that are the dynamical analogs of the Inverse Participation Ratios (I.P.R.s) in quantum localization models (2(| 

P q {L) = f d d r\^(r)\ 2q (11) 
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As a consequence of the normalization of Eq. [9l one has the identity Y q —\{L) = 1. The localization/delocalization 
transition can be characterized by the asymptotic behavior in the limit L — > oo of the Y q (L) for q > 1. In the localized 
phase T <T C these moments Y q {L) converge to finite values 

Y q (L = oo) > for T <T C (12) 

In the delocalized phase, the spreading of the polymer involves the Brownian exponent £ = 1/2, and, space being of 
dimension d — 3, the decay of the moments follows the scaling 

Y q (L) ~L-( q -V dC = L-('?- 1 )i for T > T c (13) 

Exactly at criticality, the typical decay of the Y q {L) defines a series of generalized exponents r(q) = (q — l)D(q) 



Y q typ (L) = e lny ^| T=Tc ~ L- T{q) = L- {q - 1)D{q) (14) 

The notion of multifractality corresponds to the case where D(q) depends on q, whereas monofractality corresponds 
to D(q) — est as in Eq. (|13[). The exponents D{q) represent generalized dimensions [lj : -D(O) represent the dimension 
of the support of the measure, here it is simply given by the space dimension D(0) = d = 3; -D(l) is usually called 
the information dimension [l[ , since it describes the behavior of the 'information' entropy 

s L = -^Wi(r)hi^(r) - -d q Y q (L)\ q=1 ~ Z>(l)lni (15) 

r 

Finally D(2) is called the correlation dimension [l[ and describes the decay of 



Y* VP (L) = e lnY2 ^\ T = Tc =s L _D(2) (16) 

In the multifractal formalism, the singularity spectrum /(a) is given by the Legendre transform of r(q) [1[ via the 
standard formulas 

q =/'(«) (17) 

r(q) =aq-f(a) (18) 

The physical meaning of f(a) is that the number A/l(oO of points r where the weight WL(r) scales as L~ a typically 
behaves as 

Af L (a)ocL f W (19) 
So the Legendre transform of Eq. (fTg]) corresponds to the saddle-point calculus in a of the following expression 

Y q typ (L) - J da L~ qa (20) 

The general properties of the singularity spectrum /(a) are as follows ,1| : it is positive f(a) > on an interval 
[a min , a max ] where a m i n = D(q = +oo) is the minimal singularity exponent and a max = D(q = — oo) is the maximal 
singularity exponent. It is concave /"(a) < 0. It has a single maximum at some value cxq where /(«o) — D(g = 0) 
(so here /(an) = 3), and contains the point ai = /(ai) = D(l). 

Following [H, many authors consider that the singularity spectrum has a meaning only for /(a) > fl9l I20L I23I [241. 
[25j . However, when multifractality arises in random systems, disorder-averaged values may involve other generalized 
exponents [43|, [3, IH, HE] than the typical values (see Eq. [T4l) . In quantum localization transitions, these exponents 
were denoted by f(q) = (q - l)D(q) in [H, [H 

T^(X)|t=t c =s ^" f(9) = L-^- 1 ) 6 ^ (21) 

For these disorder averaged values, the corresponding singularity spectrum f(a) may become negative f(a) < 
IH IM E B E Ell t0 describe rare events (cf Eq. [T9|) . The difference between the two generalized exponents sets 
D{q) and D(q) associated to typical and averaged values has for origin the broad distributions at criticality [2l], 
as we now describe. 
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C. Probability distributions of the Y q (L) 

The scenario proposed in (2ll . [26j in the context of quantum localization models is as follows : the probability 
distribution of the logarithm of the Inverse Participation Ratios of Eq. [TT] becomes scale invariant around its typical 
value Elm, i.e. 



\nP q {L) =]xiP q (L)+u (22) 

where u remains a random variable of order O(l) in the limit L — ► oo. According to [26[ the probability distribution 
Gl(u) generically develops an exponential tail 

Goo(u) ~ e~ x « u (23) 



As a consequence, the ratio y = P q (L) / P^ vp (L) = e u with respect to the typical value P q yp (L) — e lnPq ^ presents 
the power-law decay 

n {* < 24 > 

The conclusions of [2l|, [26| are then as follows : for q small enough q < q c , the exponent satisfies x q > 1, and the 
corresponding generalized dimensions coincide D(q) = D(q). However, for larger values q > q c , the exponent x(q) 
may become smaller x(q) < 1, and then the corresponding generalized dimensions differ D(q) ^ D(q) : the decay of 
the averaged value P q (L) is then determined by the finite-size cut-off of the power-law tail. In this case, the averaged 
values are not representative but are governed by rare events. 

In this paper, we show that this scenario for the I.P.R.'s statistics at Anderson transitions describes well our data 
for the directed polymer at criticality. 



III. RESULTS FOR THE GENERALIZED EXPONENTS D(q) AND D(q) 
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FIG. 1: (Color online) (a) Multifractality at criticality (T = 0.79) : generalized dimensions D(q) (O) and D(q) (□) associated 
to typical and disorder averaged values (Eq. 1141 and t !H) (b) Monofractality in the high-temperature phase (T = 2.): D(q) = 
D{q) = | (see Eq. [13]) 



We show on Fig. Q] (a) our results for the generalized exponents D{q) and D(q) governin g th e decay of typical 
and disorder averaged values (Eqs. [T4l and |2T|) . In agreement with the scenario proposed in [2lL [26| for Anderson 
transitions, we find that there exists a threshold q c , of order q c ~ 2 here, such that 



D{q) = D{q) for q < q c 
D{q) > D(q) for q > q c 



(25) 
(26) 
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In particular, for q — 1, the information dimension of Eg. (|15p is 

D(l) = D{1) ~ 1.5 (27) 

and corresponds to the monofractal dimension Dt>t c (<?) = 3/2 of the delocalized phase (see Eq. [T3]) . as numerically 
checked on Fig. Q] (b) . 

For q = 2, the correlation dimension D(2) defined in Eq. (fl"6|) is found to be of order 

D(2) ~ Z)(2) ~ 1.3 (28) 

For q > 3, the values for D(g) and D(q) are clearly different, in particular for q = 3 

£»(3) - 1.1 (29) 

£>(3) - 0.9 (30) 

and for q = 10 

L»(10) -0.9 (31) 

£>(10) -0.3 (32) 

The limit g — > oo will be discussed more precisely below (see Eq. I39[) 

Finally, for g < 0, we find that Y* VP {L) and Y q (L) diverge more rapidly than the power-laws of Eqs. [T4l and |2"TI i.e. 

D{q < 0) = +oo (33) 

This can be understood in the delocalized phase T > T c where it is also true, since for the free Gaussian probability 
wt=oo{'t) ~ e~( r ) l L / L?/ 2 at T = oo leads to exponential divergence of Y q in the negative region q < 0. This finding 
for our 'dynamical' model is thus very natural, but is a major difference with the Anderson localization cases where 
the generalized exponents D(q) are finite for q < 0. 



IV. RESULTS FOR THE SINGULARITY SPECTRUM f(a) 




FIG. 2: (a) Singularity spectrum /(a) : starting at some a m in = D(q = +oo) ~ 0.77 where /(a m i n ) = 0, it is tangent to 
the diagonal a — f(a) at «i = D(l) ~ 1.5 and asymptotically goes to /(+oo) = D(0) — 3. (b) Corresponding curve a(q) : 
diverging at q — ♦ 0, it goes through the point a(q = 1) = -D(l) — 1.5 and tends to amin = = +oo) as q — > oo. 

To measure the singularity spectrum /(a), we have followed the method of the q-measures proposed in [47j | . As 
explained above, in the negative region q < 0, our model does not lead to powerdaw (see Eq. I33|) . As a consequence, 
the maximum (ao, f(ao)) of the curve f(a) which is at finite distance in Anderson localization models, is rejected 
towards infinity in our case 



ao = +oo 
/(a ) = D(q = 0) = 3 



(34) 
(35) 
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Our result for the curve f(a) are shown on Fig. [2] (a) : it begins at some a m i n — D(q = +00) where f(a m in) = 0, 
it is tangent to the diagonal at ai — D(l) ~ 1.5 and asymptotically goes to /(+00) = D(Q) = 3. 

The q- values used in our computations, and the corresponding a(q) (see Eq. [20]) are shown on Fig. [2] (b). 

To measure better the minimal exponent a m i n — D{q = +00) where f{a) vanishes /(a m j„ = 0) = 0, we have 
studied the statistics of the maximal weight in each sample (Eq. [8]) 

w ™ ax =maxp[w L {f)] (36) 

The disorder-averaged value of its logarithm gives 

lnw^ ax 0.77 hxL (37) 

The first moment involves a similar value 

wZ~L-°- 75 (38) 
so our conclusion is that the minimal exponent in a typical sample is around 

a min = D(q = +00) ~ 0.77 (39) 



V. RESULTS FOR PROBABILITY DISTRIBUTIONS 



A. Probability distributions of the Y q (L) 




In Y„ 

2 u 

FIG. 3: (Color online) Histogram of lnY 2 (L) at criticality (T c = 0.79) (a) Probability distribution PL(lnY 2 ) for L = 
6,12,18,24,36,48,60,72,84,96 (b) Rescaled distributions In G L (u = \nY 2 (L) - lnY 2 (L)) : the exponential tail of Eq. (|4T|) 
is clearly visible, the corresponding slope being X2 ~ 1. 

To understand the difference between the generalized exponents associated to typical and averaged values (Eq. |2"6"| . 
we now consider the probability distributions of Y q (L) over the samples. Our results for the histograms of \yiY2(L) 
for various L at criticality are shown on Fig. [3] (a) . Remarkably, as L grows, this distribution simply shifts along the 
cc-axis with a fixed shape, as also found in [2l|, [26( for I.P.R.s at Anderson transitions. As in Eq. [22j we may therefore 
write 

\nY q (L) = \nY q (L) + u (40) 

where u is a random variable of order 0(1) in the limit L — * 00. The probability distribution G*2(w) of u = XnYziL) — 
In Y"2(£) is shown on Fig. [3] (b) for various L. ft clearly develops an exponential tail as u — » 00 

G£^oo(u) ~ e~ x - u (41) 
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FIG. 4: (Color online) Histogram of lnY2(L) in the high-temperature phase ( T = 2 ) (a) Probability distribution of Pi,(ln Y2) 
pour L = 6,12,18,24,36,48,60 (b) The distributions G L (u) of the rescaled variable u = (lnY 2 (L) - lnY 2 (L))/[AlnY 2 (L)] 
become Gaussian. 



As stressed in [2ll l26j|. the ratio y = Y q (L) / Y q vp (L) — e u then presents the power-law decay 

oc -J- (42) 
V Y* VP (L) J y 1+x * v 1 

Whenever x q < 1, the corresponding generalized dimensions differ D(q) ^ D(q) (Eq. |2!>|) : the decay of the averaged 
value Y q (L) is then determined by the finite-size cut-off of the power-law tail. Our results for the histograms of Y q (L) 
for q = 3,4, .. are similar to the results shown for q = 2 on Fig. [3] 

For comparison, we show on Fig. 0] (a) the histogram of lnY2(X) in the delocalized phase at T — 2 : as L grows, 
the width shrinks around the averaged value lnl2(i) ~ — (3/2) lnX. The corresponding rescaled distribution shown 
on Fig. |4] (b) tends towards the Gaussian distribution. 



B. Probability distributions of the last-monomer entropy sl = — w L(r) In wl{t) 




Since the last-monomer entropy is closely related to the Y q (L) (Eq. I15p . we have also computed its histogram 
over the samples both at criticality (Fig [5]) and in the delocalized phase at T = 2 > T c (Fig. [S]). Again, the rescaled 
distribution is Gaussian for T > T c (Fig. b), and strongly asymmetric at criticality (Fig. [5]b). 
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VI. FINITE-SIZE SCALING IN THE CRITICAL REGION 



For Anderson transitions, finite-size scaling involves the multifractal spectrum at criticality but a single correlation 
length exponent v (see the reviews [HI, HH). In this section, we thus try the following finite-size scaling form in the 
critical region 



Y q (L, T) = $ ( (T - T C )L X I V 



For T < T c , the convergence to finite- values Y q (L = oo, T) in the L — > oo limit yields 

Y q (L = oo,T) = (T c - TfU> with /3(g) = vf{q) 



(43) 



(44) 



This relation between the multifractal ex pon ents f(q) and the critical exponents (3{q) and v is well known for the 
Anderson transitions (see the reviews [23j, |24(). Our results for T < T c are shown on Fig. [7] for q = 2 (a) and q = 3 
(b) with the value v = 2. This value is one of the two values v ~ 2 and ~ 4 found previously in the literature for 
other observables [H, S3] . 
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FIG. 7: (Color online) Finite-size scaling (T < T c ) of = L f(q) Y q as a function of x = (T - T c )L 1/l/ , see Eq. g3] with the 
value v = 2 and for the sizes L = 12(Q), 24(D), 36(0), 48(A), 60(>) (a) q = 2 (b) q = 3 



For T > T c , the results of the finite-size scaling form of Eq. [43] are shown on Fig. [8] for q = 2 with the value v = 2. 
Note that the matching between the finite-size scaling form of Eq. 23] and the asymptotic behavior in the delocalized 
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FIG. 8: (Color online) Finite-size scaling (T > T c ) of m = L f(2) Y2 as a function of x = (T - T c )L xlv with the value v = 2 and 
for the sizes L = 12(Q), 18(*), 24(D), 36(0), 48(A), 60(t>), see Eq. gU 



phase (Eq. [T"3"]) yields a diverging amplitude -A(T) in Eq. [13] for q > 1 

y 9 (L,T>r c ) - ff^n with ~ ( T - Tc)- 1 '^- 1 )^-- 6 ^ (45) 

Our data thus points towards a correlation length exponent v ~ 2 above and below T c , i.e. towards a value very 
close to the general lower bound v >2/d = 2 of disordered systems [48j . 



VII. CONCLUSION 



In this paper, we have found that the directed polymer in a random medium of dimension 1 + 3 exhibits mul- 
tifractal properties at the critical localization/delocalization transition. We have numerically studied the statistics 
of the Y q (L) (see Eqsl8l and flO]) . which are the dynamical analogs of the Inverse Participation Ratios of Anderson 
localization quantum models [23l, [2~il . |26| . Our results are very close to the Evers-Mirlin scenario [U, [2(| for the 
Anderson transitions case. In particular, we have found that the generalized dimensions D(q) and D(q) for typical 
and disorder averaged values coincide for q < q c ~ 2 but differ for q > q c , and that the probability distributions of 
y = Y q (L)/Yq VP (L) over the samples becomes scale invariant with a power-law tail l/y 1+Xq . We have also measured 
the corresponding typical singularity spectrum f(a), which starts at the value a m i n = D(+oo) ~ 0.77, and ends at 
czmax = +oo. Off-critical results lead, through a finite size scaling analysis, to a value v ~ 2 for the correlation length 
exponent on both sides of the transition. 

Finally, our numerical results, in particular the scale invariant shape of the histogram of lnl2(^) shown on Fig. 
(a), strongly support the equality T c = T^id = 3) (see Eq. [5] and the corresponding discussion). 

Since the directed polymer in a random medium can be mapped onto a growth model in the Kardar-Parisi-Zhang 
universality class[3(|, multifractality is also expected to show up at the critical point of these growth models. More 
generally, the present study confirms that it may be interesting to characterize the critical points of quenched disor- 
dered models by their multifractal properties. 

Acknowledgements 

We thank F. Igloi for drawing our attention to many relevant references. 



[1] T.C. Halsey, M.H. Jensen, L.P. Kadanoff, I. Procaria and B. Shraiman, Phys. Rev. A 33, 1141 (1986). 
[2] G. Paladin and A. Vulpiani, Phys. Rep. 156, 147 (1987). 
[3] H.E. Stanley and P. Meakin, Nature 335, 405 (1988). 



10 



A. Aharony and J. Feder Eds, Fractals in Physics, Essays in honour of B.B. Mandelbrot, North Holland (1990). 
P. Meakin, Fractals, scaling and growth far from equilibrium, Cambridge (1998). 

D. Harte, " Multifractals, Theory and Applications", Chapman and Hall (2001). 

B. Duplantier, "Conformal Random Geometry", Les Houches, Session LXXXIII, 2005, Mathematical Statistical Physics, 
Eds A. Bovier et al., 101, Elsevier (2006). 

A.W.W. Ludwig, Nucl. Phys. B 330, 639 (1990). 

J.L. Jacobsen and J.L. Cardy, Nucl. Phys., B515, 701 (1998). 

T. Olsson and A.P. Young, Phys. Rev., B60, 3428 (1999). 

G. Palagyi, C. Chatelain, B, Berche and F. Igloi, Eur. Phys. J, B13, 357 (2000). 

C. Chatelain and B. Berche, Nucl. Phys., B572, 626 (2000). 
N. Sourlas, Europhys. Lett. 3, 1007 (1987). 

M.J. Thill and H.J. Hilhorst, J. Phys. I 6, 67 (1996). 

G. Parisi and N. Sourlas, Phys. Rev. Lett. 89, 257204 (2002). 

J. Kisker and A. P. Young Phys. Rev. B 58, 14397-14400 (1998) ; F. Igloi, R. Juhasz, and H. Rieger Phys. Rev. B 61 11552 
(2000). 

F. Wegner, Z. Phys. B 36, 209 (1980). 

C. Castellani and L. Peliti, J. Phys. A 19, L429 (1986) 

M. Schreiber and H. Grussbach, Phys. Rev. Lett. 67, 607 (1991); H. Grussbach and M. Schreiber, Phys. Rev. B 51, 663 
(1995). 

T. Terao, Phys. Rev. B 56, 975 (1997). 

A. Mildenberger, F. Evers, and A. D. Mirlin Phys. Rev. B 66, 033109 (2002) 

A. D. Mirlin, Y. V. Fyodorov, A. Mildenberger, and F. Evers Phys. Rev. Lett. 97, 046803 (2006) 
M. Janssen, Int. J. Mod. Phys. 8, 943 (1994); M. Janssen, Phys. Rep. 295, 1 (1998). 

B. Huckestein, Rev. Mod. Phys. 67, 357 (1995). 

C. Chamon, C. Mudry, and X.-G. Wen Phys. Rev. Lett. 77, 4194-4197 (1996) ; H. E. Castillo, C. Chamon, E. Fradkin, P. 
M. Goldbart, C. Mudry, Phys. Rev. B56, 10668 (1997). 

F. Evers and A.D. Mirlin, Phys. Rev. Lett. 84 , 3690 (2000) ; A.D. Mirlin and F. Evers, Phys. Rev. B 62, 7920 (2000). 
W. Pook and M. Janssen, Z. Phys. B 82, 295 (1991) 

B. Huckestein and L. Sweitzer, Phys. Rev. Lett. 72, 713 (1994); B. Huckestein and R. Klesse, Phys. Rev. B 55, R7303 

(1997); T. Brandes, B. Huckestein and L. Schweitzer, Ann. Physk 5, 633 (1996). 

F. Evers, A. Mildenberger and A.D. Mirlin, Phys. Rev. B 64, 241003 (2001). 

T. Halpin-Healy and Y.-C. Zhang, Phys. Repts., 254, 215 (1995). 

B. Derrida and H. Spohn, J. Stat. Phys., 51, 817 (1988). 

M. Mezard, J. Phys. (France), 51, 1831 (1990). 

P. Carmona and Y. Hu, Prob. Th. and Rel. Fields 124, 431 (2002); P. Carmona and Y. Hu, |math.PR/0601670| 
F. Comets, T. Shiga and N. Y oshida, Bernoulli 9, no. 4, 705 (2003). 
V. Vargas, |math.PR/0603233| 

J. Z. Imbrie and T. Spencer, j. Stat. Phys. 52, 609 (1988). 

J. Cook and B. Derrida, J. Stat. Phys. 57, 89 (1989). 

J.M. Kim, A.J. Bray and M.A. Moore, Phys. Rev. A44, R4782 (1991). 

B. Derrida and O. Golinelli, Phys. Rev. A41, 4160 (1990). 
M. R. Evans and B. Derrida, J. Stat. Phys. 69 , 427 (1992). 

C. Monthus and T. Garel, Phys. Rev. E 74, 011101 (2006) 
C. Monthus and T. Garel, Eur. Phys. J. B 53, 39-45 (2006) 

B. Mandelbrot, Physica A 163, 306 (1990) B. Mandelbrot, J. Stat. Phys. 110, 739 (2003) 
A.B. Chhabra and K.R. Sreenivasan, Phys. Rev. A 43, 1114 (1991). 
M.H. Jensen, G. Paladin and A. Vulpiani, Phys. Rev. E 50 , 4352 (1994). 

T.C. Halsey, K. Honda and B. Duplantier, J. Stat. Phys. 85, 681 (1996); T.C. Halsey, B. Duplantier and K. Honda, Phys. 
Rev. Lett. 78, 1719 (1997). 

A. Chhabra and R.V. Jensen, Phys. Rev. Lett. 62, 1327 (1989). 

J. T. Chayes, L. Chayes, D.S. Fisher and T. Spencer Phys. Rev. Lett. 57, 2999 (1986). 



